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UTILIZATION  OF  QUANTUM  DISTRIBUTION  FUNCTIONS  FOR  ULTRA-SUBMICRON 
DEVICE  TRANSPORT 

•  •• 

C.J.  Iafrate,  B.L.  Grub in  and  D.K.  Ferry 

U.S.  Army  Electronics  Technology  and  Devices  Laboratory ,  Fort  Monmouth, 
Neo-Jersey  07703,  U.S.A. 

*Scientific  Research  Associates,  Inc.  P.0.  Box  498,  Glastonbury,  Connecticut 
06033,  U.S.A. 

** Colorado  State  University,  Ft.  Collins,  Colorado  B0S23,  U.S.A. 

-  La  representation  de  Uigner  est  consider 6e  comme  une  mSthode  utilisable 
pour  decrire  les  ph£nom£nes  de  transport  quantique  dans  la  physique  des  disposi- 
tifs  submicroniques.  la  fonction  de  distribution  de  Wigner  est  rappel£e  ici  et 
sa  non  unicit£  est  discut£e,  ainsi  que  sa  generalisation  aux  autres  f one t ions  de 
distribution  quantique.  L‘ equation  du  mouvement  pour  la  fonction  de  distribution 
de  Uigner  est  dScrite  et  un  nouveau  r£*ultat,  les  moments  de  cette  equation  du 
mouvement,  est  obtenu.  On  montre  que  les  equations  du  moment  contiennent  des 
corrections  quantiques  par  rapport  aux  equations  classiques  et  que  ces  tenses 
quantiques  ne  sont  pas  negligeables  lorsque  les  longueurs  de  transit  des  porteurs 
sont  de  l'ordre  de  la  longueur  d'onde  de  de  Broglie. 

Abstract  -  The  Wigner  'representation  is  considered  as  a  method  for  describing' 
quantum  transport  phenomena  in  connection  with  submicron  device  physics.  The 
Uigner  distribution  function  is  reviewed,  its  non-uniqueness  discussed;  Its 
generalization  to  other  quantum  distribution  functions  is  examined.  The  equa¬ 
tion  of  motion  for  the  Wigner  distribution  function  is  described  and  a  new 
result,  the  moments  of  this  equation  of  motion,  is  derived.  It  is  shown  that 
the  moment  equations  contain  quantum  corrections  to  the  classical  moment  equa¬ 
tions;  these  quantum  terms  are  seen  to  be  non-negllglble  when  the  carrier 
transit  lengths  are  of  the  order  of  the  particle  deBroglle  wavelength. 

Introduction  -  As  semiconductor  technology  continues  to  pursue  the  scaling-down  of 
1C  device  dimensions  into  the  su.micron  (less  than  ten  thousand  Angstroms)  regime* 
many  novel  and  Interesting  questions  will  emerge  concerning  the  physics  of  charged 
particles  in  semiconductors.  One  of  the  more  Important  topics  to  be  considered  is 
that  of  the  appropriate  transport^  picture  to  be  used  for  a  given  spatial  and  tem¬ 
poral  regime.  Moreover,  from  the  point  of  view  of  device  physics,  it  is  most 
desirable  to  have  a  microscopic  description  of  semiconductor  transport  which  is  compu¬ 
tationally  manageable  or  at  least  amenable  to  phenomenological  treatment  so  that  its 
properties  can  be  meaningfully  incorporated  into  device  simulations.  In  this  paper 
attention  is  focused  on  a  useful  transport  methodology  for  the  ultra-submicron 
regime.  At  present,  electronic  devices  with  transit  lengths  in  the  ultra-submicron 
region  are  coming  to  fruition  due  to  the  advent-  of  MBE  processing  so  that  the  need 
for  an  appropriate  description  of  transport  in  this  regime  is  Indeed  imperative. 

Semiconductor  transport  in  the  ultra-submicron  regime  approaches  the  category 
of  quanta  transport.  This  is  suggested  by  the  fact  that  within  the  effective  mass 
approximation  the  thermal  deBroglle  wavelength  for  electrons  in  semiconductors  (See 
Fig.  1)  is  of  the  order  of  ultra-submicron  dimensions.  Whereas  classical  transport 
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physics  Is  based  on  che  concept  of  a  probability  distribution  function  which  Is 
defined  over  the  phase  space  of  the  system,  in  the  quantum  formulation  of  transport 
physics,  the  concept  of  a  phase  space  distribution  function  Is  not  possible  lnasauch 
as  the  hon-comoutatlon  of  the  position  and  aoaentua  operators  (the  Heisenberg  uncer¬ 
tainty  principle)  precludes  the  precise  specification  of  a  point  in  phase  apace. 
However,  within  the  matrix  formulation  of  quantum  mechanics.  It  is  possible  to  con¬ 
struct  a  "probability"  density  matrix  which  is  often  Interpreted  as  the  analog  of 
the  classical  distribution  function. 


There  is  yet  another  approach  to  the  formulation  of  quantum  transport  based  on 

3 

the  concept  of  the  Wigner  distribution  function  (WDF)  .  This  formalism  is  particu¬ 
larly  attractive  for  use  in  ultra-submicron  device  transport  in  that  It  contains  all 
of  the  quantum  mechanical  information  about  the  state  of  the  system  yet  has  elements 
of  the  classical  picture  implicitly  built  in.  Thus,  the  intent  of  this  study  is  to 
explore  the  potential  usefulness  of  the  WDF  as  well  as  other  possible  quantum  dis¬ 
tribution  functions  for  describing  quantum  (ultra-subalcron)  device  transport  in 
semiconductors’.  To  this  end  we  first  review  the  salient  features  of  the  WDF  and 
then  discuss  a  new  result,  the  derivation  of  the  first  three  quantum  moment  equa¬ 
tions  using  the  WDF.  It  is  shown  that  the  moment  equations  contain  quantum  correc¬ 
tions  to  the  classical  moment  equations;  these  quantum  terms  are  non-negliglble  when 
the  transit  lengths  are  of  the  order  of  the  carrier  deBroglie  wavelength. 

Wigner  Distribution  Function  -  The  Wigner  distribution  function^  is  generally  defined 
in  terms  of  all  the  generalized  coordinates  and  momenta  of  the  system  as 


pw(\ 


♦ 


pn>r  ihnrf-a, +  Jr 


(1) 


However,  for  simplicity  here  in  this  paper  we  discuss  the  properties  of  a  single 
coordinate  and  aoaentua  WDF: 
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(2) 


where  4(x)  represents  the  state  of  the  systen  In  the  coordinate  representation. 
Although  we  will  he  treating  the  WDF  for  the  special  case  of  pure  states,  the  adap- 
t  tlon  of  this  forma  lisa  to  include  mixed  states  is  accomplished  through  the  gener¬ 
alization 


pw(*.P>=  l  Pf,[  T ) 

n  —cd 


(3) 


where  P  is  the  probability  to  be  in  state  "n".  (For  example,  for  a  system  in 
contact  with  a  heat  bath  at  constant  temperature,  Pn  “  e  ^1/  .) 

The  distribution  function  of  Eq.  (2)  has  interesting  properties  in  that  inte¬ 
gration  over  all  momenta  leads  to  the  probability  density  in  real  space;  while  inte- 

3 

gration  over  all  coordinates  leads  to  the  probability  density  in  momentum  space  : 


f05  PyJx.pltfp  =  f*  Pw(*.  p)  d-k  -  ^>*(p)</>(P) 

J-a>  J-a> 


where 


<f>lp)  =  -  1  —  fm  e-|Pll/,,  \f/(x)dx 


(Aa) 


(4b) 


It  follows  from  Eqs.  (4a, b)  that,  for  an  observable  w(x,p)  which  is  either  a 
function  of  the  momentum  operator  alone,  the  position' operator  alone,  or  any  addi¬ 
tive  combination  therein,  the  expectation  value  of  the  observable  is  given  by 
<  W  >  =  jj  w  Pw(x,  p)tfx<^p  which  is  analogous  to  the  classical  expression  for  the 
average  value.  Herein  lies  the  interesting  aspect  of  the  WDF;  it  is  possible  to 
transfer  many  of  the  results  of  classical  transport  theory  into  quantum  transport 
theory  by  simply  replacing  the  classical  distribution  by  the  Wigner  distribution 
function.  However,  unlike  the  density  matrix,  the  WDF  itself  cannot  be  viewed  as 
the  quantum  analog  of  the  classical  distribution  function  since  it  is  generally  non- 
positive  definite  and  non-unique. 

Further  resemblance  of  the  Wigner  distribution  function  to  the  classical  distri¬ 
bution  function  is  apparent  by  examining  the  equation  of  time  evolution  for  Pw(x,p). 
Upon  assuming  that  *(x)  in  Eq.  (2)  satisfies  the  Schrodinger  equation  for  a  system 
with  Hamiltonian  H”p^/2m+V(x) ,  it  has  been  shown  that  P^  satisfies  the  equation 


dPw  p_ 
dt  m 


+  B  ■  Pw  «  0 


(5) 
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where 

J_  ®  » (li/2)2"4,1  aa"»‘v(K)  3MlPjl*,p) 

6t**~  n  »,ol"  (2n+l)!  axi*+i  dp*"*'  <«) 

Alternately,  d.P^  cen  be  expressed  ae  ~  { "jjj"  V{*)Pw(*,p) 

where  1C  la  understood  chat  Che  position  gradient  operates  only  on  the  potential 

energy,  V(x).  In  the  limit  1rM3,  Q  .  p  .  =  -  -4^  so  that  Eq.  (6)  reduces  to  the 

w  dx  dp 

classical  collisionless  Boltzmann  equation. 

The  Wigner  distribution  function  is  not  the  quantum  analog  of  the  classical 
4 

distribution  function  as  evidenced  by  the  fact  that  it  is  generally  non-positive 
defnite  and  non-unique.  The  non-uniqueness  is  clearly  evident  by  noting  that  many 
quantum  distribution  functions  can  be  constructed  which  obey  the  sum  rules  of 

4  5 

Eqs.  (4a, b).  One  other  distribution  function  which  we  have  examined  in  detail  is  * 


Pw  - 


2wft 


f^dy  ^4  *-y>e,py/K  s 
-co 


•/lirtT 


**(x)eip*£(p> 


(7) 


where  4(p)  is  given  by  Eq.  (4b).  In  actuality  the  non-uniqueness  and  non-positive 
definiteness  of  the  Vigner  distribution  function  can  be  traced  to  the  non-commutation 
of  the  position  and  momentum  operators. 

Moment  Equations  -  In  this  section,  we  derive  the  first  three  moments  of  a  "tfigner- 
Boltzmann'-llke  transport  equation, 


_P_  j£w 
dt  ™  i* 


(8) 


port  equation,  since  P^  is  not  a  true  probability  distribution  function. 


This  equation  was  constructed  to  include  an  ad-hoc  collision  term  which  stay  not 
necessarily  express  the  same  phenomenology  as  that  of  the  classical  Boltzmann  trans- 

These 

problems  are  conceptually  reduced  when  dealing  with  moments  in  a  relaxation  approxi¬ 
mation.  The  moment  equations  are  obtained  by  multiplying  Eq.  (6)  by  an  appro- 
priate  function  of  momentum,  d(p),  and  then  integrating  over  all  momenta  to  obtain: 
d  <  <£>  _  I  d 

~rrT 

(9) 

. 

where  «  >  refers  to  an  Integration  over  momentum.  In  making  the  assumptions  that 
4(p)  he  an  analytic  function  of  momentum  and  that  P  (x,p)  vanish  at  the  momentum 
limits,  it  follows  that 


dt 


dx 


i  d2"*'  ,  r ~  ** 

«?o\  2i  /  (2n+0!  dx*"+l  V(*U  ^p)  $p 


f*{ p)  (l^r  pw)^p  /(  ?<p>)  pw*p 


(10) 


and  that  Eq.  (9)  becoaes 


I  4  r? 

- -  ♦  ~=r  -  < d>p> 

41  m  4x  Y 


9  (jlT  -J _ ( vu)y.^  >  =  <274*) 

■£UJ  (2»+lH \4*2»+l  7  4p2,"M  *\  dt  4.11. 


(11) 


For  specific  values  of  T(p),p01Pi  end  Eq.  (11)  becomes 

+  JL  -L<P>  =0. 

41  in  4x 


4<i*  + 
41 

1 

m 

4 

“3T 

<p2  >  +  p 

4V  . 
4x 

<  P  1 

_L  A. 

2  - 

i 

4  , 

<P> 

4v 

pm  41 

<P  > 

+ 

2m2 

- -  <  P  ^ 

dx 

+  m 

dx 

'  COll. 


where  upon  reduction 


(12a) 


(12b) 


<P  >  \2')  f0('°  (n-j)!J!  dx1  4x"'i 

thereby  showing  the  dependence  of  <p°> .  For  n“0,l,2,3  we  show  <pn>  explicitly: 


(13) 


<p°>  =  y/*  f 


-*£*) 


(14s) 


<  p2  > 


'  \  2 1/  \  ^  2x*  2  4x  4x  4x  V  (14b) 

<d3>  .  (*$(+*£*  -3  *L*  (140 

<P  >  -  ^  ax5  3  dx  4x2  4x2  4*  4x5  V 


Note  that  use  of  Eqs.  (3Ad  in  Eq.  (L2a)  results  in  the  correct  quantua  mechanical 
continuity  equation. 

In  order  to  see  the  transition  between  the  quantua  and  classical  regimes,  we 
invoke  a  wavefunction  of  the  fora  ^(x,1)  *  A(x,1)e  .and  so. 
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where  p(x,t)  le  Che  probability  density  and  v(x,t)  if  the  ensemble  velocity.  Insert¬ 
ing  the  vavefunctlon  of  Eq.  (15)  into  Eq.  (14)  results  in 

<p°>  *  <p*»  *  (16a) 

<p*>  *  (mv)*/)  -  ln/>  (16b) 

<  ps>  =  ( mv)J/>  -  p  {  3mv  \n/>  *  (mv)}  (lk) 

It  i«  evident  fro*  Eq.  (16)  that  the  tens  possessing  sn  explicit  dependence  on 
"B"  are  the  quantum  corrections  to  the  classical  nonen tun-density  soseots.  Further¬ 
more,  insertion  of  these  nonentun-noaents  Into  Eqs*  (12)  results  in  a  set  of  ■osmt 
equations  which  contain  explicit  quantua  corrections  as  veil*  Of  course9  as  tl"*0» 
these  movent  equations  reduce  to  the  classical ,  sero  temperature,  sosent  equations. 

The  logarlthalc  derivative  tern  appearing  in  Eqs.  (16)  is  not  unconnon;  such  a 
tern  appears  in  the  real  part  of  the  Schrodlnger  equation  when  the  vavefunctlon 
associated  with  Eq.  (15)  is  used6.  In  the  limit  where  this  term  is  negligible,  the 
Schrodlnger  equation  goes  over  to  the  Bamilton-Jacobl  equation. 

We  have  made  order  of  magnitude  estimates  of  the  strengths  of  the  quantum  correc¬ 
tion  tern  appearing  in  Eq.  (16).  Using  a  Gaussian  spatial  variation  for  P(x),  and 
a  thermal  ensemble  valuefor  ymv*,  typical  of  central  valley  GaAs  electrons,  we  find 
that  the  quantum  correction  is  substantial  for  widths  of  the  order  of  SO  A. 
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